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GCN: a gaseous Galactic halo stream? 
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ABSTRACT 

We show that a string of Hi clouds that form part of the high-velocity cloud complex 
known as GCN is a probable gaseous stream extending over more than 50° in the Galactic 
halo. The radial velocity gradient along the stream is used to deduce transverse velocities as 
a function of distance, enabling a family of orbits to be computed. We find that a direction of 
motion towards the Galactic disk coupled with a mid-stream distance of ~ 20 kpc provides 
a good match to the observed sky positions and radial velocities of the Hi clouds comprising 
the stream. With an estimated mass of 10 s M Q , its progenitor is likely to be a dwarf galaxy. 
However, no stellar counterpart has been found amongst the currently known Galactic dwarf 
spheroidal galaxies or stellar streams and the exact origin of the stream is therefore currently 
unknown. 
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1 INTRODUCTION 

The study of neutral hydrogen (Hi) structures in the Milky Way 
dates back to the early 1950s, with the first detections of 21cm 
emission by Ewen & Purcell (1951) having followed the predic- 
tion by van der Hulst that this line should be detectable astro- 
nomically. Since these initial days of radio astronomy, the Galactic 
Hi sky has been mapped repeatedly, providing very detailed maps 
of the Hi distribution in the Milky Way as projected on the sky 
(Leiden-Dwinge loo Survey: [ Hart mann & Burton IT997I . Hi Parkes 
All Sky Survey: Barnes et alj2001 , Institute Argentino de Radioas- 



tron omfa Survey: Bajaja et al 



vey: 



2005, Southern Galactic Plane Si 



McClure-Griffifh s et all 120051 , Arecibo Galactic Hi Survey: 



IStanimirovic et all 120061) . However, detections of gaseous struc- 
tures in themselves do not, unfortunately, provide us with infor- 
mation about their three-dimensional spatial distribution. In the 
disk plane, it is possible to de-project the sky-projected distribution 
from knowledge of their radial velocities and by assuming rotation 
around the Milky Way (e.g. for the Galactic spiral arm structure; 
iLevine et al]|2006h . In most cases, however, pinning down the lo- 
cation of a gaseous structure requires the knowledge of a stellar 
population that it can be associated with. 

High- velocity clouds (HVCs) of neutral hydrogen are a class of 
Hi structures whose kinematic behaviour cannot be reconciled with 
that expected for disk-plane material. While several theories exist 
for their origin (the Galactic fountain, strip ping from dwarf galax- 
ies, remnants of galaxy formation; see e.g. Bregman 2004), none 
of these have outright popularity over others, primarily due to the 
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fact that HVCs must be treated on a case-by-case basis. When de- 
ducible, their distances and metallicities can help to constrain their 
origin, but the former is particularly difficult to obtain for many of 
the Galactic HVCs; their locations at often high Galactic latitudes 
means that identifying and utilising a halo field star to apply the in- 
terstellar absorption line technique, from which distance limits can 
be derived, is usually a taxing task. 

In the context of stripped gas associated with dwarf galaxies, 
the Sagittarius dwarf galaxy is the only Milky Way satellite, aside 
from the Magellanic Clouds, for which there is an Hi feature that 
could have been s tripped from the main body of the dwarf galaxy 
dPutman et aHl2004l) . Indeed, the HVCs comprising the Magellanic 
stream are the only ones for which there is definitive knowledge 
that the gas originated from the Milky Way satellite companions 
themselves. In terms of satellite galaxies in the Local Group, the 
Magellanic system is unique in the sense that the two galaxies 
have clearly been interacting with one another as well as with the 
Milky Way; the binary interaction has most probably been crucial 
in producing the gaseous stream of over 100° in length. There 
are two reasons why we might expect to see more examples of Hi 
streams than this one system that is currently known. The first is 
that studies derived from cosmological simulations indicate that 
a significant fraction of dark matter subhalos f all into their more 
massive host halo in groups <Li & Helmill2008h . In this scenario, 
dwarf galaxies that are hosted within such subhalos will naturally 
have companions that may be involved in close interactions, 
depending on the degree of binding within the infalling group. The 
second reason for expecting halo Hi streams is that dwarf galaxies 
venturing too close to the parent galaxy should have their gas 
content at least partially stripped through tidal forces and/or ram 
pressure. Indeed, the Hi content of Milky Way and M3 1 satellites 
are well known to correlate with distance from their respective 
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parent galaxies jGrcevich & Putmarl2 009). 



In this letter, we highlight the possibility that a collection of 
clouds belonging to an HVC complex known as GCN may be part of 
yet another gaseous stream at an average heliocentric distance of ~ 
20 kpc. In Section[2] we first outline how the velocity of the stream 
can be calculated as a function of distance. Section [5] presents the 
orbit derived from this analysis and the results are summarised in 
Section!] 



2 DERIVING THE TRANSVERSE VELOCITY OF A 
STREAM 

Before turning to the analysis of the Hi stream, we recall the tech- 
nique for deducing the transverse velocity at a given location along 
a stream by using t he gradient in line-of-sight velo city there, as ha s 
been described by I Jin & Lvnden-Belll d2007h and lBinnevI d2008h - 
Its strength lies in the use of line-of-sight velocity data along an 
extended structure, so that proper motion measurements are not re- 
quired to constrain an orbit, if at least one distance is known. If no 
distance information is available, then the transverse velocities are 
calculated as distance-dependent quantities that can then be used to 
derive a family of possible orbits in a given gravitational potential. 
Depending on the morphology of the orbit, both the distance of the 
stream and its direction of motion can be constrained by comparing 
the resulting orbit with observational data. 

The formalism for deducing the component of the velocity 
vector transverse to the line of sight at a given point on a stream 



was provided in a previous paper jjin & Lvnden- Bell 2007). Here, 
we simply quote the solution to the resulting quadratic equation for 
this transverse velocity, v 3 : 



(dvi/dx) -4(l.V<!/> )d 



(1) 



where vi is the line-of-sight velocity corrected to the Galactic Stan- 
dard of Rest (GSR), 1 is the normalised line-of-sight vector, d 
is the heliocentric distance, ip is the Galactic potential and \ is 
the angle along the stream, measured from some fiducial point. 
The full velocity vector at any point along a stream is given by 
v = vi\ + VsS, where s is the unit vector in the plane of the sky, 
along the apparent direction of motion of the stream at 1; one pos- 
sible method by which s may b e calculated is given by equation 
(16) in Ijin & Lvnden-Belll d2008h . The two solutions for v s usually 
provide orbits in opposite senses. Here, a comparison of the com- 
puted orbits with the Hi data turns out to be sufficient to specify the 
direction of motion, given the large angular extent of the gaseous 
structure considered. 

In the following analysi s, we us e a th ree-component 
Galactic model described by Pacz vnskil Jl990l) . which com- 
bines a Miyamoto- Nagai model for the disk and spheroid 
(Mivamoto & Na gail[l975h . and a near-logarithmic, spherical po- 
tential for the halo. We also briefly describe the effect of changing 
the halo potential in Section [3Tl 



3 RESULTS 

3.1 Determining the orbit of GCN 

Using the all-sky HVC map by IWakkeJ J2004 figure lb) as 
a visual guide, we extract several GCN Hi features in the 



X (deg) 

Figure 1. Variation in line-of-sight velocity along the HVC complex GCN, 
where \ is 'he angle along the stream as measured from the point of lowest 
latitude and increases towards the Galactic plane. The error bars on the ve- 
locities indicate the width of the velocity peak in the LAB data. The dashed 
line indicates the velocity gradient at the starting point for the orbital inte- 
gration. 






Figure 2. Orbital evolution in Galactic longitude and latitude (£, b) and 
line-of-sight velocity relative to the Galactic Standard of Rest (u;). The 
panels provide a comparison of GCN data (dots) with an orbit (solid line) 
computed using parameters derived from its properties at I = 39.3°, 
b = —31.0° (see main text for details). The computed orbit places the 
Hi complex at 25 kpc at the IC point. The open circle represents the CHVC 
at b = 47° (see discussion in Section (3/2). Note that this solid-line orbit 
does not take this point into account; its inclusion here serves merely to il- 
lustrate the position of the CHVC relative to the stream's orbit. The dotted 
line describes an orbit that is constrained to have the same IC point as the 
orbit shown by the solid line, but passes closer to the CHVC in I, b and vi . 
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Figure 4. Distance evolution of the derived orbit for GCN (corresponding to the orbit indicated by solid lines in Figures[2]and[3j, shown in Hammer projection. 
The open squares denote positions of the Hi clouds comprising GCN, while the open circle represents a compact high-velocity cloud that may or may not be 
associated with the main stream. Triangles indicate the locations of Milky Way dwarf spheroidal galaxies, as listed in Table[T] 



angle along stream (deg) 

Figure 3. Variation of heliocentric distance with angle along the stream 
(measured from the IC point), for the orbits shown in Figure [2] The solid 
line describes the case in which the CHVC is not included and corresponds 
to the orbit shown in Figure [4] while the dotted line is with its inclusion. 
The central region of the solid line, in black, corresponds to the extent of 
GCN that we have followed in the LAB data. 



Leiden/ Argentine/Bo nn (LAB) Galactic Hi survey data cube 
jKalberlaet all [2005) lying along the general direction of I ~ 
33 — 40° with negative velocities relative to the gsrQ Figure Q] 
shows the evolution of vi along the string of HVCs that indicate a 
strong communal trend in both sky position and velocity, the data 
for which can then be used to calculate dvi / dx and hence the trans- 
verse velocity along the stream. 

In defining the initial condition (IC) point for the orbit calcu- 



1 Another HVC complex (GCP) overlaps GCN spatially for a good chunk 
of the latter. However, these two structures have rather contrasting line-of- 
sight velocities so that any physical association between them is clearly 
ruled out and their emission is easily separated in the LAB data. 



lation, it is important to choose a location on the stream where the 
behaviour of dvi/dx is well defined. This usually corresponds to 
a point near the middle of the stream's extent. Here, the IC point 
chosen to satisfy these criteria is at a Galactic longitude and lati- 
tude of (£, b) — (39.3, —31.0)°, with the corresponding properties 
Vi — — 147.5kms~ 1 saddvi/dx = —3.0 kms _1 /deg (indicated 
by the dashed line in Figure [TJ. The free parameter in this analy- 
sis is the heliocentric distance to the stream. However, there is also 
some uncertainty in the velocity gradient along the stream of order 
0.5 kms _1 /deg; we therefore treat this gradient as a somewhat 
variable parameter. In order to recover a good match to the position 
and velocity of the Hi data points, we find that the distance to the 
IC point is reasonably well constrained at ~ 25 kpc with the re- 
maining orbital parameters as stated above0 The orbit that results 
from taking this set of parameters is one of the best found and is 
shown by the solid lines in Figure[2] where the three panels indicate 
the variation of each of the variables I, b and vi as a function of the 
other variables. Note that the stream mainly extends in Galactic lat- 
itude (in other words, x is effectively given by b), and so the panel 
showing the trend of Vi in longitude is less instructive in charac- 
terising and understanding the orbit than the other two panels. The 
solid line in Figure [3] shows the evolution of heliocentric distance 
along this orbit, with the highlighted central section corresponding 
to the extent traceable in the LAB data. 

We note that the details of any orbit calculated are natu- 
rally dependent on the choice of the Galactic potential. As many, 
m ore cosmologically -motivated halo profiles than that described 
by IPaczvnskil dl99fj|) now exist in the literature, the analysis was 
also re-performed with Paczyh ski's halo being replaced by an 
adiabatically-contracted NFW dNavarro et al. 1996, 1997) profile 



2 Altering this initial distance by ±10 kpc leads to orbits that are clearly 
discrepant with one half of the stream. 
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Table 1. Properties of Milky Way dwarf spheroidal galaxies shown in Fig- 
ure|4] The columns give the Galactic longitude and latitude, heliocentric ra- 
dial velocity, line-of-sight velocity relative to the Galactic Standard of Rest 
and heliocentric distance. Radial velocities and distances are taken from 
those compiled by iMateol Jl998t> and lMartin et alj J2008|). with radial ve- 
loci ties for the 'ultra-faint ' dwarf galaxies taken from lMartin et alj J2007I) 
and Simon &Gehdi2007l) . with additional informatio n for Segue 1 , Leo V , 
Seg ue 2 and Pisces II fromlGeha et alj j2009l) . lBelokurov et al jl2008Ll2009l) 
and lBelokurov et alj f20irjK respectively. 



with parameters constrained by |Xueetal.ld2008h . in order to check 
the dependency of the results on the choice of the halo potential. 
Keeping the same initial conditions as before, we find that the 
change in halo potential leads only to small changes in the result- 
ing orbit that are mainly manifest above the Galactic plane. Within 
the region of sky covered by GCN, the differences in the results 
are negligible in both the radial velocity trend and sky projection 
of the orbit. The distances are slightly affected by the fact that the 
Xue et al. halo contains less mass than Paczyhski's. As a result, the 
orbit has a larger apocentre and reaches a smaller pericentre, with 
a maximum change in distance of ~ 10% compared to the origi- 
nal orbit. The qualitative results and most of the quantitative results 
are, however, unchanged by the modification of the halo potential. 



3.2 Checking for associations 

A different representation of the stream's orbit is provided in Fig- 
ure[4] where the 'anti-centre projection' allows us to directly com- 
pare our results with the Wakker (2004) HVC map. The same orbit 
as that shown by the solid lines in Figures [2] and [3] is now overlaid 
with the Hi stream data points as open squares; additionally, all of 
the known Milky Way dwarf spheroidal (dSph) galaxies are plotted 
as filled triangles (Table [TJ. Although there are a few spatial coin- 
cidences, the distance and velocity of the orbit at these locations 



firmly rule out possible associations of this Hi stream with any of 

the currently known Milky Way dSphs. 

The Milky Way is also host to ~ 150 globular clusters dHarrisI 

1996). Of these, a handful of globular clusters with latitudes be- 
tween —50° and 0° lie close to the stream but, with velocity dif- 
ferences of more than 100 km s -1 from the GCN clouds or the 
inferred orbit, connections to the Hi stream are unlikely. The same 
conclusion holds for three other globular clusters at 40° < b < 50° 
that lie within 10° of the orbit there, although one should remember 
that the orbital path at these latitudes is already a strong extrapola- 
tion of the orbit derived for the main part of the Hi stream residing 
at negative latitudes. 

Thanks to the advent of large sky surveys, numerous detec- 
tions of Galactic stellar streams have been reported in the literature 
in the last few years. We therefore also considered the possibil- 
ity that the GCN stream might have a counterpart in one of these 
recent discoveries, if not in the globular clusters or intact dwarf 
spheroi dals. We find that no known stellar str eam [Monoceros 
stream dNewberg et all |2002|; llbata et alj [2003b . Orphan stream 
I Grillmaij2006l:lBelokurov et alj2007h7 Grillmair-Dionatos stream 
l Grillmair & Dionatosll2006l : IWillett et ai]|2009l) ; Acheron, Cocy- 
tos, Lethe, and Styx streams (Grillmair 2009); solar neighbour- 
hood st reams in SPSS DR7 dKlement et alj l2009h; Cetus polar 
stream dNewberg et alj [2009b ; Pisces overdensity dWatkins et alj 
2009)] coincides with the orbit determined in this letter for the 
GCN stream. Although inhabiting a similar region of the sky, the 
stream is also offset from the Magellanic system in location and 
velocity. Combining this with our check for coincidences with the 
stellar streams therefore leads us to conclude that the Hi stream has 
no currently known stellar counterpart. 

Finally, we find that this orbit passes very close to a compact 
high- velocity cloud (CHVC) in the Wakker HVC map at (£, b) = 
(18, 47)° . Perusal of the LAB data shows this CHVC to have a line- 
of-sight velocity of (—107 ± 15) kms" 1 . As shown by the dotted 
lines in Figures [2] and [3] its inclusion in our analysis helps to con- 
strain the distance of the orbit further by requiring the orbit to trace 
as much as possible the clouds over a larger extent on the sky, but 
the association of this Hi clump with the rest of the structure is un- 
clear from the current data. Higher resolution Hi dat^fl would be 
required to see whether this CHVC exhibits an elongation along the 
direction of motion expected from the orbit; naturally, detection of 
lower column density material along this direction would also aid 
in solidifying this hypothesis. 

3.3 Discussion 

Despite its 'orphan' nature in having no obvious progenitor or stel- 
lar counterpart, our findings provide evidence for the first halo Hi 
stream outside of the Magellanic stream at distances where it is 
conceivable for Hi structures to display relatively ballistic motion. 
Although the orbits were always integrated for several wraps, here 
we have shown only the first 660 Myr to either side of the IC point. 
The reason for this is three fold: the first is simply a matter of clar- 
ity, so that the reader does not have to untangle an array of different 
orbital wraps on such a map. The second is that the exact form of 
the gravitational potential is, as always, somewhat uncertain. In ad- 
dition, we do not expect the motion of any halo stream (whether 

3 For example, ver y recent results from the Effelsberg-Bonn Hi survey 
i Winkel et al. 2010) clearly show head-tail features and filamentary struc- 
tures in new, high resolution data for GCN clouds. 
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stellar or gaseous) to follow an exact orbit, as the extended pro- 
file itself is created through changes in orbital energy and angular 
momentum amongst the stream's constituents. This is even more 
important for gaseous structures than it is for stellar streams, as gas 
is subject to non-gravitational influences in the form of ram pres- 
sure and drag forces, effects that do not concern a stellar stream. We 
therefore do not expect to find a perfect orbital match and, rather, 
content ourselves with an orbit that reproduces well the general 
trend in position and velocity of the string of Hi clouds of our in- 
terest. 

As we were only able to identify a few of the brightest re- 
gions in the LAB data, it is difficult to make an accurate mass 
estimate of the stream. We the refore simply use a representative 
column density of 10 19 cm -2 JWakkeill2004l) and a total angular 
size of 20 deg 2 at an average distance of 20 kpc to arrive at a con- 
servative lower mass estimate of 2 x 10 5 Mp. Studies of the H i 
content of Local Group dwarf galaxies dGrcevich & Putmar]|2009l) 
indicate this to be a plausible gas content for an undisturbed, dis- 
tant dwarf galaxy. Observing such a quantity of Hi spread into a 
stream would, however, clearly be more challenging than observ- 
ing it as a bound and more concentrated structure, and may explain 
why the stream had not been identified before. It is also not sur- 
prising that the GCN stream would be much less massive than the 
Magellanic stream given the masses of its progenitors, the Large 
and Small Magellanic Clouds, which themselves are a few orders 
of magnitude more massive than most of the dwarf companions of 
the Milky Way. 



4 SUMMARY 

We have shown that a string of Hi clouds that form part of the 
high- velocity cloud complex known as GCN is a probable gaseous 
stream in the Galactic halo at a distance of ~ 20 kpc. We deter- 
mine its orbit by utilising the large gradient in line-of-sight velocity 
along its extent, and determine its direction of motion to be towards 
the Galactic plane. We also identify a compact HVC, whose sky po- 
sition nearly intersects with the forward projection of the stream's 
orbit, and whose velocity differs by approximately 50 km s" 1 from 
that of the orbit there. This association, however, remains a spec- 
ulative one before further Hi studies along the direction of the or- 
bit presented here can be performed. The origin of such a gaseous 
stream with an estimated gas content of a few times 10 5 M© is most 
probably a satellite galaxy of the Milky Way. However, no progen- 
itor has been found amongst the currently known Galactic dSphs 
or globular clusters through a comparison of their locations and ve- 
locities with the stream's orbit, hence the exact origin of the stream 
is currently unknown. 
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